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Abstract 

We calculate the cross section for the production of Ty-mesons via 
AiV -^ NNrj in a relativistic One-Boson-Exchange-Model. Using this 
cross section we then determine the probability for the production of 

m . 

f^ ' an Ty-meson by a A-resonance moving in nuclear matter. The result is 

f^ ' compared to prescriptions in BUU-calculations in which 77-production 



proceeds both through a direct channel and through the sequential 
process A -^ Nn; ttN — > Nrj. 
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1 Introduction 

The detection of mesons produced in heavy-ion-collisions is an important 
means for the investigation of the behavior of nuclear matter under extreme 
conditions. An analysis of pion spectra, for example, has shown that within 
the reaction zone, where at bombarding energies of about 1 GeV per nucleon 
baryon densities up to 3po are reached, about a third of the nucleons is 
excited to A-resonances 0]. Numerical simulations of heavy-ion-collisions 
arrive at the same result [^. 

A direct consequence of this resonance matter formation is in these cal- 
culations an enhancement of sub-threshold particle production P, § . If the 
bombarding energy per nucleon is below the free threshold for the produc- 
tion of a particle, even when Fermi motion is taken into account, then a 
large part of the observed reaction cross section is, according to numerical 
simulations, due to entrance channels involving resonances; these resonances 
store energy which can then be used to produce this particle. 

While the calculated particle yields generally agree quite well with exper- 
iment when the resonance channels are included, it is clear that this result 
depends on the cross sections used for the production processes involving 
resonances in the entrance channel; the BUU-calculations mentioned just 
use an ad-hoc recipe for the direct reaction channel AA^ — > AAmeson and 
contain an explicit treatment of the sequential channel A -^ Ntt; ttN — > Nij. 

It is thus the aim of this study to investigate the elementary production 
process in more detail in order to put the results on the dominance of res- 
onance channels on a firmer ground. For this study we chose the process 
AA^ -^ NNt] because all the interaction vertices occurring there are known 



from other processes. 

2 The process AN -^ NNr] 

We use a relativistic One-Boson-Exchange-Model to calculate the cross sec- 
tion of the process AA^ -^ NNr] . Models like this have been used success- 
fully to describe the production of rj-uiesons in collisions between nucleons 

[5-7]. 

2.1 Couplings 

We use the standard interaction-terms in the hadronic Lagrangian. Because 
the A has isospin 3/2, the only exchange bosons that have to be taken into 
account are the isovector mesons vr and p. 

The interaction-terms involving the pion are (the subscript N* always 
refers to the A^(1535)-resonance): 

CANn = ^^^ ^'^T^Nd^n +h.c. 

T 
^NN'Tv = S'AfAf'TT ^iV* 2*^ ^ + /i.e. . (1) 

For the interactions containing the p-meson we use: 

Canp = i ^^''' ^^ f Y 75 *JV {d^P. - d,p^) + h.c. 
niA + riiN 

Cnn'p = igNN*p ^Af*757M2^ P^ + ^■'^- ■ (^) 

Finally the coupling of the rj is taken to be: 

CNN*-n = gNN*r) ^N*'^N rj + h.c. . (3) 



The diagrams that are to be considered are shown in Fig. 1 (Top). We ne- 
glect pre-emission graphs since these are suppressed kinematically. Graphs, 
where the rj is not produced via a A^(1535) but via a nucleon, are neglected as 
well, because the NNrj-veitex is assumed to be negligible due to an analysis 
of Bhalearo and Liu |8|. In principle, there could also be a A^(1535)-A- 
vr-coupling. However, according to [^ the branching ratio for the decay 
A^(1535) -^ An is smaller than 1%. Therefore this coupling is clearly negli- 
gible. 

For the off-shell particles in the diagrams in Fig. 1 (Top), we use 
monopole form factors: 

/(,^) - ^ (4, 

2.2 Parameters 

The value of the A-nucleon-pion coupling is easily determined using the 
decay-width of the A: 

QAN-K = 2.15 (5) 

which corresponds to a A-width of 115 MeV. This number is rather fixed 
because of the well known width of the A. 

The experimental values of the A^(1535)-width and its branching ratios 
have big uncertainties P; in order to fix gNN^T^ and gNN*r) we therefore 
consider the process vrA^ — > N-q. These coupling-constants were already 
determined by other groups (for example HP]). However, different models 
or different energy ranges were employed in these analyses so that we are 
forced, in order to ensure consistency within the model used, to fix these 
numbers using our parameterization. 



To calculate the cross section for vrA^ -^ Ntj one has to evaluate the 
diagrams shown in Fig. 1 (Bottom). The couplings gNN*rj and gNN*TT enter 
in two different ways into these diagrams: First, they appear as factors 
at the vertices and, second, they enter through the A^(1535)-width in the 
propagator. The dependence of the A^(1535)-width and its branchingratio 
on these couplings was determined by evaluating the diagrams for the decay 
into a pion and an rj, respectively. The two-pion decay of the A^(1535) was 
assumed to give a 10 % contribution to the total width. 

Evaluating the diagrams in Fig. 1 (Bottom) numerically and comparing 
these calculations with the data available for ttN -^ Nr] we arrived at the 
numbers (Fig. 2): 

9nn*tt = 1-37 

gNN*r, = 1-53 (6) 

The couplings involving the p-meson are harder to determine. However, the 
following calculations show, that the contribution of the p-exchange is much 
less important than the exchange of a pion. We, therefore, took the numbers 
of other groups (e.g. Q) for these couplings: 

QANp = 13.3 
gNN'p = 0.61 (7) 



(see Sec. 12^). 

2.3 Kinematics, the exchanged pion 

The kinematics of the reaction AA^ -^ NNij need further attention, because 
the exchanged pion can be on-shell. This can be seen from the diagrams 



in Fig. 1 (Top): If the exchanged meson is a pion and one simply cuts 
the meson line, one ends up with two diagrams that are both possible as 
reactions between on-shell particles. The first one stands for the A-decay 
and the second one for the process vrA^ -^ Nrj. A detailed analysis shows, 
that the pion can be on-shell for cm-energies of the AA^-system from about 
50 MeV above the threshold for r/-production up to about 7 GeV. In this 
case the pion propagator becomes infinite, and the integration over the phase 
space of the final particles cannot be performed because the integral diverges. 

The solution to this problem is to remember that the reaction considered 
here does not take place in the vacuum, but during a heavy-ion-collision, i. 
e. within a medium with a non-vanishing baryon density. In such a medium 
the pion acquires an in-medium width which renders the pion propagator 
finite, so that the total cross section can be calculated. 

This width enters into the propagator via the pion selfenergy S^: 

The relation between the width and the imaginary part of the pion selfenergy 
is: 

ImS^ = -S^r, . (9) 

The width of the pion is nothing else but its inverse lifetime and can therefore 
be expressed by means of the mean free path: 

T^ = - = -^ = v^ pa^N + ^nip) , (10) 

r ( 

where Vn is the pion velocity, p is the nuclear density and anN is the to- 
tal pion-nucleon cross section, which contains elastic scattering as well as 



inelastic processes [10|. By using the form ([To|) we assume a local-density 



approximation for the selfenergy of the pion; we are aware of the fact, that 
this is a good approximation only in the nuclear bulk matter region. 

r,r (p) in Eq. ( p!o| ) represents the part of the selfenergy of the pion which is 
due to many-body effects in the nuclear medium. The effect of 2- and 3-body 
absorption on the pion selfenergy in the region of the A-resonance has been 
discussed extensively in the literature (e.g. [12, lH|). In the present case how- 



ever the pion-nucleon system has an invariant mass above the //-threshold 
where several resonances contribute (N(1440), N(1520), N(1535),. . . ), so 
that there is no quantitative information available. In the absence of any 
reliable information on the many-body part of the pion selfenergy we have 
therefore simply multiplied the vrA^ cross section with a density dependent 
factor by making the replacement 

CTnN > 0-7rAr(l + a{p)) . (11) 

with a{p = 0) = 0. Since we will work at constant density (local density 
approximation), effectively this just amounts to multiplying the cross section 
(TttAt with a constant factor. 

The expression for the imaginary part of the selfenergy in Eq. (^ can 
be derived by evaluating the lowest order diagrams contributing to Stt by 
means of Cutkosky's rule. For a = it reduces to the well known form 
Im Stt = — wivr Ttt if one chooses the rest frame of the pion as the reference 
frame. In our case the more general form (^ must be used, because the 
width given by ( [l0| ) uses the rest frame of the medium as the reference 
frame. 

For the real part of the selfenergy we took results of calculations em- 
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ploying the A-hole-model [pi] ]. The A-hole- model also yields an imaginary 
part of S which is similar to ( p^ ) |jl^. However, a.„j\f then only contains 
the A and no other resonances. We improved this by using a parameteriza- 
tion of the experimental data for a-j^N in (P^, thereby including also higher 
resonances. 

An important consequence of Eq. ( p!o|) is that because of the density 
dependence of the pion-width the cross section for AA^ -^ NNr] depends on 
the density as well. This will be discussed further in the next section. 

2.4 Results 

The results of our calculations for the total cross section are given in Fig. 
3 (Top) for three different nuclear densities and a of Eq. ( pT| ) set to zero. 
Interpreting the density dependence of the cross section at fixed ^/s one 
has to bear in mind, that a cross section in general is a transition amplitude 
normalized with respect to the flux of the incoming particles and the number 
of target particles. So the cross section for A~p — > nnr] as plotted in Fig. 3 
does not contain the probability for a A meeting a nucleon, which does, of 
course, also depend on the density. 

Keeping this in mind the density dependence of the cross section at fixed 
^/s is easy to understand. It becomes smaller with increasing p, because the 
probability for the exchanged pion to react with the medium grows with 
the density. In the -y/s-region where the pion can be on-shell (see Sec. |2.3| ) 
the total cross section scales roughly with 1/p and diverges for vanishing 
density. 

Fig. 3 (Bottom) shows the cross section with and without the exchange 
of a p-meson. One sees that the process AA^ -^ NNr] is clearly dominated 



by the exchange of the pion (The couphng constants used were given in Eq. 
(0)). A more recent analysis |l^ finds a A-A^-/9-couphng which is about a 
factor of two below the one given in Eq. (0), so that the contribution of the 
p-exchange is even smaller. The large contribution of the /o-exchange found 
for A^A^ -^ NNr] in [^, ^ does not appear here, because the on-shell pion 
dominates the reaction. 

3 Comparison to the BUU-Parameterization 

The purpose of this work is to provide microscopically justified input for 
transport theories like the BUU-model to describe the production of r]- 



mesons in heavy ion collisions. In standard BUU-calculations |14] the pro- 
cess AA^ -^ NNt] is implemented in two additive ways: 

1. Sequential production: A A-resonance can decay into a nucleon and 

a pion, that is treated as a real, on-shell particle. One of the possible 
reactions for this pion is then vrA^ -^ Nij. 

2. Direct production: One of the reactions allowed for the A besides the 

decay is the direct production of an rj-uieson via AA^ -^ NNr] . The 
cross section for this process was assumed to be equal to the cross 
section for NN -^ NNr] at the same ^/s, corrected for the different 
spin-isospin degeneracies in the entrance channel. 

In the last section it was shown, that the total cross section is governed 
by the exchange of an in-medium pion with a finite width. Consequently the 
contribution to the r/-production by on- and off-shell pions cannot strictly be 
separated. Indeed, Salcedo et al. [p!5|] have shown for a similar problem that 



in the limit of vanishing density the sequential and the direct production 
processes give exactly the same result. Thus, a BUU-calculation employing 
both the sequential and the direct production channel might contain some 
double counting. This is the point we want to investigate in the following 
discussions. 

In order to do so we consider a single A-resonance of given energy moving 
in homogeneous matter with density p and calculate the r/ production rate 
by solving appropriate rate equations. For simplicity we take one definite 
isospin channel, namely a A^ moving within a medium of protons. In this 
channel the additional cross sections we will need below are experimentally 
well known. Also for simplicity we first set Tt^{p) of Eq. ( p^ equal to 
zero, thus neglecting many-body reactions of the pion; these effects will be 
discussed separately at the end of this section. 

The rate equations have to contain the decay widths and reaction rates 
of all possible processes. These are for the 

A channels: 

A-decay: This enters into the rate equations via the decay width Fa. 
In our scenario there is no medium modification of this width due 
to Pauli blocking, because a A^ can only decay into a vr^ and 
a neutron, for which there is no Pauli blocking in a medium of 
protons. 



A-absorption: Given the cross section for AA'^ -^ NN [14] the reac- 
tion rate for this reads: 

FAp^nn = "UA ^Ap^nn P (12) 

where va is the velocity of the A and p is the density of the 
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protons. 
Direct r/-production: The reaction rate for this is given by: 

^Ap^nnrj = Vj\ CT/^p^nnrj P ■ (13) 

So far, no sequential production is included. To account for it we have 
to consider the reactions that the pion produced by the decay of the A can 
undergo: 

Secondary vr interactions: 

Reaction with the medium: Using the total pion-nucleon cross sec- 
tion, which appeared already in Eq. (p!o|), the reaction rate for 
this is: 

r,rp = ■Vtt CTttp P (14) 

As has been said before, being the total cross section a-j^p contains 
elastic and inelastic processes. 

ry-production: Our calculation for the cross section for up — > nrj are 
shown in Fig. 2 together with experimental data. With that cross 
section the reaction rate is: 

Combining all these channels we arrive at the following rate equations: 

— N/\{t) = (-TAp^nn - Ta - FAp^nrHj) Nl^,{t) 

— iV^(t) = rAp^„n,,iVA(t)+Frp^n,,A^7r(i) • (16) 
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These equations constitute a first collision model; as such they do not con- 
tain any refeeding of the A or of the vr. This is not only sufficient for our 
purpose of comparing the BUU-prescription with our microscopic calcula- 
tion, but is even a cleaner test because it does not sum over complicated 
collision histories with their inherent averaging over many processes and 
cross sections. 



The first of Eqs. (IC) contains some double counting, because on-shell 
pions appear both in Ta and TAp^nnrj- However, for normal nuclear den- 
sity it is TAp~*nnri ~ 0.01 Fa, so that this double counting has negligible 
consequences. 

Because of the absence of any refeeding terms in the A and vr channels 
Eqs. ([iq) can easily be solved. If one uses the initial conditions iVA(O) = 
1; Ntj{0) = Nfj{0) = and takes the limit P^ = linii^oo ^riit), then P^ is the 
probability for the production of an r/-meson in the first reaction of either 
the A or the pion produced in the decay of that A with a proton of the 
medium. 

If the density dependent cross section calculated above is used in Eq. 
(|l3|) for T^p^nnrj, the width r^p^„^ has to be set equal to zero, because 
then the cross section (JAp^nnr) already contains all kinematically possible 
invariant masses of the intermediate pion. The results of this calculation 
are shown in Fig. 4 (upper curves). The density dependence of P^ is weak 
because the approximate 1/p-dependence of the cross section is multiplied 
with the p-dependence of the collision rate (Eq. p!3|) . 

Also shown in Fig. 4 is the result for the sequential process, obtained 
by setting (TAp^nnr? to zero in the rate equations given above (lower curves). 
For p = pq the sequential r/-production probability is about 30 % smaller 
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than the direct production cross section given by the upper curves. For 
decreasing density, however, the upper and the lower curves converge and 
for p ^ they coincide, in agreement with the theoretical results of Salcedo 
et al |l^. The physical reason for this is that for /? ^ the pion has to 
travel an arbitrarily long distance before reacting with a nucleon. This can 
only happen, if the pion is strictly on shell, so that the sequential description 
becomes exact, and the upper and lower curves in Fig. 4 have to coincide. 
The difference between the two sets of curves thus reflects an important 
in-medium effect. 

It can also be seen from Fig. 4, that at finite density the direct re- 
production and the sequential process have different thresholds. This is 
due to the fact that the sequential threshold lies above the absolute energy 
threshold for rj production in the AA^ channel because the requirement of an 
on-shell pion constitutes an additional constraint on the reaction dynamics 



(see Sec. 2^). Only for p = both thresholds agree, in agreement with the 
theoretical considerations of Salcedo et al. 

The Fermi momentum of the protons is not taken into account in Fig. 4. 
Including this leads to a smearing of all curves over a certain energy range, 
which blurs the differences shown in Fig. 4. 

In order to check the validity of the parameterization of the direct pro- 
duction cross section and the inclusion of both direct an sequential pro- 
duction channels in standard BUU-calculations we now compare our results 
with calculations that treat the possible reaction channels like the BUU- 
model. In the context of our scenario this means that we have to include in 
the rate equations both the sequential production as in Fig. 4 and the direct 
production process, but now with the density-independent parameterization 
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given in |14] for (7Ap~*nnr}- This does not change the result in the hniit of 
p — > because in this case T/\p^nnr} vanishes while Fa remains constant 



(first of Eqs. 16). The two contributions together with the resulting total 
production probability are shown in Fig. 5. 

Fig. 6 (Top) shows a comparison of the results of the BUU-simulation 
including both the direct BUU-type production and the sequential channel 
(i.e. of the results of Fig. 5) with the results of the calculations employ- 
ing our density dependent cross section for the direct channel only. The 
two curves are remarkably similar although there is some difference in their 
energy dependence. This difference is mainly due to a somewhat different 
energy dependence of the calculated cross section for AA^ — > NNt] and its 
BUU-parameterization. Taking the Fermi motion into account diminishes 
this difference significantly (Fig. 6 (Bottom)). The agreement between the 
two prescriptions at other densities than po is similar. 

So far, all the calculations have neglected the contribution of many-body 
effects on the pion selfenergy. Fig. 7 now shows results for two different val- 



ues of a{po) of Eq. ([Tl]), entering via the pion selfenergy and Eq. (14), 
respectively. In both cases the agreement between the BUU prescription 
and the calculation with our calculated cross section is still very reason- 
able for low to intermediate kinetic energies of the A (Ta < 1 GeV) . This 
agreement for different values of a is due to the fact that the replacement 



(11) affects the results of both models in an analogous way by causing a 
higher reaction probability of the pion. Only at higher energies we obtain a 
stronger dependence of Pr^ on a{po) than in the BUU-case. This is simply 
due to the fact that in the BUU-case the cross section for the direct process 
is taken to be density-independent, so that in this case only the sequential 



14 



production depends on a. 

4 Conclusion 

We have calculated the cross section for AA^ -^ NNr] in a microscopical 
model using empirical interaction vertices. Because of the special kinematics 
of this reaction, the exchanged pion can be on-shell so that a straightforward 
evaluation of the corresponding Feynman graphs leads to divergent results. 
By including in-medium effects, however, we have obtained a finite, density 
dependent result. This result contains all kinematically possible situations, 
including the point, where the invariant mass of the exchanged pion is equal 
to its free mass. 

We have then used this cross section to check the ad-hoc treatment of 
the AA^ channel for particle production in transport theories like the BUU- 
model. This is essential, because at energies far below threshold ry-mesons 
are created mainly via the resonance channel. 

The BUU-prescription allows both the direct and the sequential process. 
The sequential production involves a pion which is propagated as an on- 
shell particle, while the direct process is included on an ad-hoc basis with a 
guessed cross section to account for the ry-production via off-shell pions. In 
order to compare our 'exact' calculation with this BUU-recipe we have set 
up a simple scenario in which the probability for the initial production of 
an ry-meson by a A of given energy was calculated for both cases by solving 
rate equations. 

As a result we have found, that by allowing sequential production alone 
one obtains an r/-yield, which was clearly below the one obtained with the 
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cross section calculated in this work (Fig. 4). Only if we included sequential 
and direct production at the same time, using the cross sections as given 
in ||l^, the correct ?7-production probability was reproduced quite well (Fig. 
6, Top). The different energy dependence still visible there vanished nearly 
completely when the Fermi distribution of the medium was taken into ac- 
count (Fig. 6, Bottom). We thus conclude, that there is no double counting 
in the BUU-treatment and that the BUU-recipe for the AA^ -^ NN-q cross 
section works surprisingly well. This result puts the enhancement of far- 
subthreshold particle production obtained in BUU-calculations on a firmer 
basis. 
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Figure 1: Top: Diagrams dominating the process AA^ — > NNij . Bottom: 
Diagrams contributing to the process vrA^ — > Nr]. 

Figure 2: Experimental data for the cross section for ir^p -^ nij together 
with our calculations like described in the text. The data were taken from 



m- 



Figure 3: Total cross section for A p ^ nnrj. Top: Assuming three different 
nuclear densities. Bottom: With and without the p-exchange (p = po). 

Figure 4: ?7-production probability for different densities using only the di- 
rect process with the density dependent cross section (upper curves) and 
using only the sequential production (lower curves) without Fermi motion. 

Figure 5: Total ry-production probability in a BUU scenario (solid curve). 
The lowest (dotted) curve gives the direct production probability obtained 
using the BUU-parameterization. The dashed curve shows the sequential 
contribution (same as in Fig. 4). 

Figure 6: Top: r/-production probability using only the direct process with 
the density dependent cross section calculated in this paper and using the 
sum of the sequential and the direct process in the BUU-parameterization 
(for p = pq). Top: Without Fermi motion. Bottom: With Fermi motion. 

Figure 7: r/-production probability for two different values of a(po) of Eq. 
(O) (with Fermi motion, p = po) 
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